#if !defined(SEQUENCE_GAPS_H)
#define SEQUENCE_GAPS_H

//888888888888888888888888888888888888888888888888888888888888888888888888888888888888888
//
//	File		 :	sequence_gaps.h
//	Contents	 :	Header file for functions handling gaps in protein sequences
//	Creation Date:	14/04/2003
//	Last Modified:	14/04/2003	lg
//	tabs width	 :	every 4 characters
//
//888888888888888888888888888888888888888888888888888888888888888888888888888888888888888

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000
#include <vector>
#include <string>



struct gap_locus
{
	unsigned beg;
	unsigned len;
	gap_locus(unsigned b, unsigned l): beg(b), len(l) {}
};

typedef std::vector<gap_locus> t_gap_loci;
bool expand_gaps(const t_gap_loci& gaps, std::string& prot);
bool expand_gaps_cdna(const t_gap_loci& gaps, std::string& cdna);
void remove_gaps(t_gap_loci& gaps, std::string& prot);


#endif
